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The magnetic properties of the weak itinerant ferromagnet ZrZn2 are analyzed using Landau 
theory based on a comparison of density functional calculations and experimental data as a function 
of field and pressure. We find that the magnetic properties are strongly affected by the nearby 
quantum critical point, even at zero pressure; LDA calculations neglecting quantum critical spin 
fluctuations overestimate the magnetization by a factor of approximately three. Using renormalized 
Landau theory, we extract pressure dependence of the fluctuation amplitude. It appears that a 
simple scaling based on the fluctuation-dissipation theorem provides a good description of this 
pressure dependence. 



The physics of metals near ferromagnetic quantum 
critical points (QCPs) has attracted renewed interest fol- 
lowing several recent discoveries of materials with un- 
usual still poorly understood transport and thermody- 
namic properties, as well as unusual low-temperature 
states, particularly superconductivity coexisting with 
ferromagnetism. 1-3 From a generic, qualitative point 
of view, the phenomena are understood as being con- 
nected with renormalization, scattering and pairing due 
to strong fluctuations in the ferromagnetic order param- 
eter (i.e., spin fluctuations) as the critical point is ap- 
proached. Besides the mentioned superconducting tran- 
sition, in clean samples of many weak itinerant ferro- 
magnetic metals the magnetic transition near the QCP 
crosses over from the second order to the weakly first or- 
der, although whether this happens in ZrZn2 is not yet 
established. The physics of this cross-over is not clear 
yet. In any case, quantitative, material specific under- 
standing of these phenomena is still lacking. 

ZrZn2 is a prototypical example of a weak itinerant 
(Stoner) ferromagnet. Very small magnetic moments 
(0.12 to 0.23 hb) have been reported. These do not 
saturate even at fields up to 35 T, indicating softness 
of the magnetic moment amplitude and suggesting exis- 
tence of soft longitudinal spin fluctuations. The Curie 
temperature, Tc drops approximately linearly with pres- 
sure, starting at « 29 K at P = and decreasing to 
ss 4 K at P — 16 kbar, 4 which extrapolates to a QCP 
at P = 18 — 20 kbar. The discovery of superconductiv- 
ity in the ferromagnetic (FM) phase 3 resulted in renewed 
interest in this compound, including several theoretical 
studies (Refs. 5-9 and others). The relative structural 
simplicity of this compound and the availability of high 
quality experimental data as functions of H, T and P on 
clean samples suggest this material as a test case for de- 
veloping understanding of quantum critical phenomena 
in ferromagnetic metals. Here we focus on the magnetic 
properties, in particular the renormalization of LDA re- 
sults due to fluctuations. 

Density functional theory is in principle an exact 
ground state theory. It should, therefore, correctly de- 



scribe the spin density of magnetic systems. This is usu- 
ally the case in actual state of the art density functional 
calculations. However, common approximations to the 
exact density functional theory, such as the LDA, may 
miss important physics and indeed fail to describe some 
materials. A well know example is in strongly Hubbard 
correlated systems, where the LDA treats the correla- 
tions in an orbitally averaged mean field way and un- 
derestimates the tendency towards magnetism. Overes- 
timates of magnetic tendencies, especially in the LDA, 
are considerably less common, the exceptions being ma- 
terials near magnetic QCPs; here the error comes from 
neglect of low energy quantum spin fluctuations. Indeed, 
the LDA is parameterized based on the uniform elec- 
tron gas at densities typical for atoms and solids. How- 
ever, the uniform electron gas at these densities is stiff 
against magnetic degrees of freedom and far from mag- 
netic QCP's. Thus, although the LDA is exact for the 
uniform electron gas, and therefore does include all fluc- 
tuation effects in the uniform electron gas, its descrip- 
tion of magnetic ground states in solids and molecules 
is mean field like. This leads to problems such as the 
incorrect description of singlet states in molecules with 
magnetic ions as well as errors in solids when spin fluc- 
tuation effects beyond the mean field are important. In 
solids near a QCP, the result is an overestimate of the 
magnetic moments and tendency toward magnetism (i.e. 
misplacement of the position of the critical point) due to 
neglect of the quantum critical fluctuations. 10,11 Exam- 
ples include Sc 3 In, 12 ZrZn 2 , 6 and Sr 3 Ru 2 7 (Ref. 13). 
Sr3Rii207 displays a novel metamagnetic quantum criti- 
cal point, 14 while, as mentioned, ZrZn 2 shows coexistence 
of ferromagnetism and superconductivity. The effects of 
such quantum fluctuations can be described on a phe- 
nomenalogical level using a Ginzburg-Landau theory, in 
which the magnetic properties defined by the LDA fixed 
spin moment curve are renormalized by averaging with 
an assumed (usually Gaussian) function describing the 
beyond LDA critical fluctuations. 15,16 Although a quan- 
titative theory allowing extraction of this function from 
first principles calculations has yet to be established, one 
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can make an estimate based on the LDA fixed spin mo- 
ment curves as compared with experiment. 

LDA calculations of the magnetic properties of ZrZn2 
(Refs. 5-7, 17) are found to be sensitive to shape approx- 
imations, possibly because of the very small energy scales 
involved. In particular, atomic sphere approximations re- 
sult in smaller magnetizations than those found by more 
accurate full potential methods. 5-7 In fact, it was found 
that full potential calculations produce a Stoner factor 
of w 1.16, as opposed to 1.01 in the atomic sphere calcu- 
lations, indicating a stronger tendency to magnetism in 
the full-potential calculations. 17 

Our full potential LDA calculations for ZrZn2 at its 
experimental volume yield a magnetization of 0.72 \ib 
per formula unit - 3 to 4 times larger than experiment, 
reflecting an crucial role for the renormalization of the 
magnetization due to beyond LDA fluctuations, presum- 
ably associated with the QCP. 

As mentioned, one can incorporate such fluctuations 
into LDA calculations by renormalizing the Landau ex- 
pansion for the free energy with Gaussian spin fluctua- 
tions of a given r.m.s. amplitude (see Refs. 15, 16 and ref- 
erences therein). The latter can be obtained empirically, 
or estimated from the parameters of the band structure 
with an ansatz to separate spin fluctuations included in 
the LDA from those neglected, as discussed in Ref. 16. 
Here we report renormalized Landau functional calcula- 
tions where one parameter, the r.m.s. amplitude of the 
beyond LDA fluctuations at P = is taken as an ad- 
justable parameter, determined by comparison with the 
experimental P = magnetization, and use it to de- 
scribe, without further empirical parameters, the pres- 
sure and field dependence of the magnetic properties of 
ZrZn 2 . 

The LDA calculations were done using the gen- 
eral potential linearized augmented plancwave (LAPW) 
method. Local orbital extensions were included to accu- 
rately treat high lying core states and avoid lineariza- 
tion errors. 18 ' 19 The Hedin-Lundqvist exchange corre- 
lation function was used with von-Barth-Hcdin spin 
scaling. 20,21 The valence and Zr semicore p states were 
treated in a scalar rclativistic approximation, while the 
core states were treated fully relativistically. LAPW 
sphere radii i?=2.1 a were employed with a dimension- 
less basis set cutoff RK max —9. Brillouin zone samplings 
were done using the special k-points method, with 182 
points in the irreducible 1/48 wedge of the zone. Con- 
vergence tests were done, showing that these parameters 
were adequate. For example, fixed spin moment calcu- 
lations at the experimental lattice parameter were done 
using up to 1300 points in the wedge, with very slight 
changes of less than 0.01 /is in the magnetization. Cal- 
culations at the experimental lattice parameter were also 
done with a different sphere radius, R— 2.45 ao, again 
with negligible changes in magnetization. 

In order to construct the Landau expansion, we did 
fixed spin moment calculations, determining the total en- 
ergy as a function of magnetic moment an volume, using 



seven lattice parameters from 13.0 ao to 13.9 oo plus the 
experimental lattice parameter of 13.9358 ao. 22 The vari- 
ation of the energy with volume yields a bulk modulus, 
B=l.O Mbar, which we use to set the pressure scale, 23 
since there is no experimental value in the literature to 
our knowledge. Using this value, the QCP at P = 18 — 20 
corresponds to a volume compression of 1.7-1.9%. 

We now turn to the magnetic properties in the LDA. 
As shown in Fig.l, the magnetization drops slowly from 
0.72 fiB at zero pressure (V=338 a 3 ,) to 0.68 [Ib at 
V= 299 a 3 ,, P = 161 kbar (all volumes and magneti- 
zations are given per formula unit). At this pressure 
the ground state becomes nonmagnetic and the moment 
suddenly collapses to zero. The ferromagnetic state re- 
mains metastable until y=290 a 3 ,, P =212 kbar. Thus, 
the LDA predicts not a QCP, but a first order transi- 
tion at a pressure of P « 161 kbar. Leaving aside the 
question of the order of the transition, the LDA strongly 
overestimates the magnetization, and has a much higher 
transition pressure than experiment, implying an overes- 
timate of Tc as well. Additionally, the LDA yields very 
weak P dependence of the moment up to the transition 
pressure, while experiment finds moments that decrease 
considerably with P until at least P m 16 kbar. 24 

To proceed, we use the fluctuation-renormalized Lan- 
dau theory. 25 A large literature exists on this subject, 
for instance, the review of Ref. 15. The basis of this the- 
ory is that the main omission in LDA calculations is from 
long-range ferromagnetic spin fluctuations, which are im- 
portant near a QCP. One writes the Landau expansion 
of the LDA total energy as 

E LDA (M) = a + J2^a 2n M 2n , (1) 

n>l 

and then introduces additional Gaussian zero-point fluc- 
tuations of an r.m.s. magnitude £ for each of the d com- 
ponents of the magnetic moment (for a 3D material d = 
3). After averaging over these, one obtains a fluctuation- 
corrected functional. The general expression 15 ' 16 reads 

E re normalized{M) = a + J2 n >l 2~^ 2 ™^ 2 "' ^ 

i>o a 

Two approaches are, in principle, possible at this point: 
one is to evaluate £ using the fluctuation-dissipation the- 
orem, 

however, this requires some knowledge of the susceptibil- 
ity x(q, u)), and a choice for the cutoff in the integration. 
This choice gives the separation between the fluctuations 
accounted for in LDA from those missing. In the most 
pessimistic view it converts one unknown parameter, £, 
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into another, though it should be said that the cutoff 
may be much less material and pressure dependent than 
£ itself. The other approach is to treat £ as an adjustable 
parameter. Here we are interested in the magnetic phase 
diagram of ZrZn 2 in a pressure range corresponding to 
that where magnetism is observed experimentally, so it 
is possible to adjust £ to reproduce the magnetic moment 
at ambient pressure and then use it for the whole pressure 
range. The fluctuation-dissipation theorem, though not 
used directly, is used implicitly to construct an ansatz 
for the P dependence of £ 2 : the lowest-order expansion 
of the bare susceptibility xoi'ii^), 

Xo (q, w) - N(E F ) - aq 2 + ibuj/q (4) 

gives rise, near a QCP, to the formula (see, e.g., Refs. 26, 
27), 

g2 = 647V(^) 2 [Q4ln(1 + Q _ 4)+ln(1 + Q4)] 



where Q = q c y/a/bvp, q c is a cutoff in the momentum 
space, is the Brillouin zone volume, and vp and N(E F ) 
are the Fermi velocity and the density of states, respec- 
tively. The expression in the square brackets depends 
on its argument logarithmically, so the main volume de- 
pendence comes from the prefactor. Following the ar- 
guments of Rcf. 16, this prefactor scales with V as V^ 1 
in the effective mass approximation. Thus, in the first 
approximation we write £ 2 (U) = £ 2 {Vo)Vo/V. 

In order to ensure stable fits, we have chosen the min- 
imal power in the expansion 1, n = 6. Fig. 2 shows the 
quality of the fits, which is quite good. The value of £(Vb) 
that yields the experimental value of the magnetic mag- 
netization, M = 0.17 /is, is then found to be £(Vb) = 0.5 
fiB- Although the resulting dependence of £ on V is rel- 
atively weak, its effect on the phase diagram is large: In 
Fig. 3 we show the (zero-temperature) equilibrium mag- 
netization in zero field, as a function of volume. One can 
see that neglecting the volume dependence of £ leads to 
a QCP at P c fts 29 kbar, while using the above scaling, 
one gets nearly exact value P c « 15 kbar. We should 
recall, however, that this is the idealized phase diagram 
in zero field, while actual measurements are performed in 
a small, but finite field. Near a QCP even a small field 
can change magnetization drastically, as Fig. 4 illustrates. 
It is interesting to note that the metamagnetism present 
in the bare LDA calculations disappears when the renor- 
malization is included account and as a result a QCP 
appears. In reality, it may be that symmetry breakings 
other than uniform ferromagnctism occur near the QCP 
and change the transition to first order. It would be very 
interesting to experimentally investigate whether this in 
fact is the case, and if so how close to the transition it 
occurs and what the relevant order parameter is. 

In summary, we report LDA calculations of the mag- 
netic energy of ZrZn 2 under pressure. Our results demon- 
strate that the LDA substantially overestimates the ten- 
dency to magnetism in the whole experimentally studied 



pressure range. This is an indication of strong quan- 
tum spin fluctuations, associated with the QCP. Us- 
ing fluctuation-renormalized Landau theory, we find that 
spin fluctuations with an r.m.s. amplitude of 0.5 [Ib are 
needed at P~0 to obtain agreement with the experimen- 
tal magnetization. We further find that using a simple 
scaling based on the fluctuation-dissipation theorem we 
are able to describe the phase diagram up to the critical 
pressure with a good accuracy. 
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FIG. 1. Unrenormalized LDA magnetic phase diagram: the 
solid line is the calculated magnetic moment for those pres- 
sures where a magnetic solution exists (left axis); the dashed 
line is the magnetic stabilization energy, (right axis, same 
scale as left axis but units are mRy) Note a metamagnetic 
behavior at P > 161 kbar: there exists a magnetic solution, 
although its energy is higher than that of the nonmagnetic 
state. 




FIG. 2. Fixed spin moment calculations for lattice param- 
eters 13.15, 13.30, 13.45, 13.60, 13.75, and 13.90 a . Solid 
lines are the 6-th power fits according to Eq. 1. 
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FIG. 3. Magnetization as a function of pressure, calculated 
from Eq.2, using either a constant £ = 0.5 [j>b or with £ 2 scaled 
as the inverse cell volume (the right curve). Dots show the 
experimental magnetization at zero pressure and the experi- 
mental critical pressure. 
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FIG. 4. Magnetization as a function of pressure, calculated 
with the scaled £ 2 , in an external fields of 0, 1, 2, 3, or 4 T. 
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FIG. 5. Magnetization as a function of field; the pressures 
are to 20 kbar, spaced by 2 kbar, with alternating light and 
heavy lines. 
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